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Abstract 

We construct a model of non-uniform condensate having a spatially 
modulated complex order parameter that makes it kinematically an x- 
ray solid, i.e., a real mass density wave, but one admitting an associated 
superfluid flow. Intrinsic to this state is a non-classical translational inertia 
which we derive for the case of a potential flow. Connection to the 
non-classical rotational inertia observed in recent experiments on solid 
helium-4 is discussed. Our semi-phenomenological treatment suggests a 
flow-induced supersolid-to-superfluid transition. 

Recent experiments of Kim and Chan [1,2] on solid helium-4 at very low tem- 
peratures have strikingly revealed a non-classical rotational inertia that seems 
intrinsic to it, much as it is in the case for superfluid He II. Such a supersolid was 
indeed predicted much earlier on theoretical grounds as a plausible concomitant 
of a quantum crystal with delocalized defects, or of a Bose-Einstein condensate 
[3-5], and had motivated years of research [6]. The non-classical inertial effect 
was, however, estimated to be very small, and direct tests were therefore sug- 
gested [5]. Thus, the question "can a solid be superfluid?" , raised some 35 years 
ago [5] has now been finally answered in the affirmative. In this work we ex- 
plicitly construct a complex order-parameter that underlies a real mass density 
wave modeling the supersolid, and demonstrate analytically its necessarily non- 
classical inertia — the signature of a supersolid. We also predict a flow-induced 
supersolid-to-superfluid transition. 
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In order to motivate an order-parameter approach to the supersolid, let us 
first recall that geometrically a solid is a periodic spatial modulation of matter 
density. Such a crystalline solid structure will be revealed kinematically in its 
characteristic x-ray diffraction peaks — we may call it an x-ray solid [3]. This 
kinematic description, however must be supplemented by the energetics of its 
stability against deformation. We will show below that both these conditions are 
realized in our order-parameter description of the supersolid admitting a non-zero 
superfluid flow relative to the laboratory frame in which the density wave is at 
rest. We will specifically consider a potential flow and show that the translational 
inertia of the supersolid is smaller than its literal mass. The potential flow is 
consistent with the geometry (topology) of the experiments [1,2] in which the 
flow is confined to the narrow annular space between two co-axial cylinders, and 
has quantized Onsager-Feynman circulation [= nh/m] taken around the annulus. 
The situation can be more complicated for a simply connected topology with 
quantized vortices in it. We will then comment on its relation to the diminished 
rotational inertia observed in the above cited experiments. 

Our construction of the model supersolid state for a system of interacting 
Bose particles such as ^He is informed by the following physical considerations. 
The interparticle interaction (repulsion) is known to deplete the condensate. The 
interaction, however, favours large fluctuations, i.e., an instability, towards long- 
range diagonal order corresponding to a crystalline solid, as evidenced by the 
peak in the liquid helium static structure factor at the roton minimum wavevec- 
tor. Failure to form a solid (under its saturation vapour pressure) is, however. 
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prevented by the large zero-point kinetic energy due to the high value of the 
de Boer quantum parameter (= h/a{mey^'^ , where a and e are, respectively, 
the range and the depth of the interaction potential). It solidifies only under 
pressure that offsets the zero-point pressure through the PV term in the Gibbs 
potential. The large zero-point amplitude comparable to the interparticle sepa- 
ration, however, persists and the associated delocalization leads to the possibility 
of an off-diagonal long-range order co-extensive with the diagonal crystalline or- 
der, i.e., the supersolid phase. With this picture in mind, consider a system of 
scalar Bose particles at zero temperature, and let its number density (n(x)) be 
modulated in space as a real density wave 

n(x) = no + ^ Ug cosg • x, (1) 
g 

where summation over g spans the reciprocal lattice vectors for the periodically 
modulated condensate. 

The complex order parameter 4>{x) underlying the above real density wave 
modulation, and obeying in general the Gross- Pitaevskii equation [7], is then 

V'(x) = ^n(x)exp(i0(x)), (2) 

where ip{x.) is as usual the macroscopic wavefunction obtained by Bose condens- 
ing a macroscopic number (N) of the Bose particles into a single one-particle 
state. Our simple model assumes complete condensation at zero temperature. 
Now, associated with this complex order parameter ^(x) there is a particle- 
number current density j(x), given by 

j(x) = * (x)VV^(x) - C.C.) = ^n(x)V0, (3) 
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In the laboratory frame in which the real density modulation is assumed to be 
pinned still. Thus, we have the stationary flow with V . j(x) = 0, which de- 
termines j(x) for the given boundary condition. The phase 0(x) can then be 
obtained by solving 

V0=(?) . (4) 

Our semi-phenomenological treatment of the supersolid phase, and of the 
possible supersolid-superfluid transition is now in principle as follows. With an 
appropriate choice for the dominant reciprocal lattice vectors for the density wave 
in Eq. (1), corresponding to the peak in the liquid-helium static structure factor 
[8], and for a given flow j(x), we calculate the kinetic energy associated with the 
underlying complex order parameter V'(x). The use of a single dominant set of 
wavevectors, of course, means considerable delocalization of the atoms about the 
nominal lattice sites. This is, however, all the more justified for the solid helium 
than for a classical solid for reasons of large zero-point amplitude for the "^He 
atoms. The part of the kinetic energy involving the flow j(x) quadratically can 
now be expressed in terms of the total mechanical momentum associated with 
the flow. The real density modulation, of course, remains at rest in the laboratory 
frame. This at once identifies the inertia associated with the flow that will turn 
out to have a non-classical value. It will involve the modulation amplitude r]g to 
be determined by the overall minimization of the Gibbs free energy. The part 
of the kinetic energy not involving the flow is to be identified as the zero-point 
energy due to the extent of localization implied by the density modulation. It will 
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contribute an important part to the Gibbs free energy, and is known to prevent 
solidification of ^He under its saturation vapour pressure. It has to be off-set by 
an external pressure. Thus, we first consider the usual free energy where we follow 
the conventional mean-field theory of liquid-solid phase transition, but add to it 
the kinetic energy associated with the complex order parameter that underlies the 
real mass density wave as discussed above. Minimization with respect to rjg then 
determines the density-wave amplitudes and the phase transition involved. The 
above procedure is straightforward in principle, but algebraically cumbersome in 
3 dimensions. The ideas relevant to the supersolid phase are, however, contained 
essentially in the 1-dimensional model in the mean-field sense. This is the case 
we will now treat analytically. 

Specializing Eqs. (1-4) to the ID case, and retaining the single dominant 
density wave for the wavevector of magnitude g [8], we have from Eq. (4) for 
the phase 0(x) [9] 



where we have set j(x) = jo = constant for the uniform stationary flow. 

The kinetic energy associated with the complex order parameter il>{x) is 
readily calculated to be 
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where E^aow is associated with the flow (jo), and can be re-written as 

p2 

-^flow = /_ (7) 



with P — Lmjo — the total momentum associated with the flow, and M — 
mnoL — the total literal mass of the system. Thus, we can identify M^l — r)l 
with the associated inertia! mass of the supersolid which is non-classical for 

M„„„_e,assical = M{1 - rilf' < M (8) 

which depends on the depth of modulation rjg (to be determined from free-energy 
minimization). Equation (8) shows that the non-classical inertia associated with 
the flow decreases monotonically with rjg from M (for rfg = 0, no solid-like density 
modulation) to zero (for rjg — 1, a solid). Such a trend is clearly reasonable 
physically. The term -Bzero-pomt in the expression for the kinetic energy Ek in 
Eq. (6) is the zero-point energy due to the extent (~ 2TT/g) of localization 
associated with the spatial modulation at wavevector magnitude g. 

Now, we turn to the free energy Fo[rjg] (actually energy at zero temperature 
here with the entropy term omitted) whose minimization should give the den- 
sity modulation rjg and thus determine the phases and the phase transition. As 
discussed above, we have in the single dominant density-wave meanfield approx- 
imation [8] 

Fo/N = roVg + Uo4vt + ^oevt + " • • > (») 
where the subscript '0' denotes no flow (jo = 0)). (Note the omission of the 
cubic term on the R.H.S. of Eq. (8), well known in the context of classical 



liquid-solid transition in the 3D case where for a dominant |g| it is possible to 
choose three vectors forming an equilateral triangle for, e.g., an f.c.c. lattice 
[8]. For a ID case, however, this is not possible). To this now we must add 
Ek for a non-zero flow as obtained above. Also, a PV term is to be added, 
where P is the external pressure the system is subjected to. This term can be 
readily shown to give —SrioP in the Gibbs free energy for our ID model, where 
S is the reduction in the lattice spacing due to P. (For the ID case here, P is, 
of course, an external force). Thus, we have for the Gibbs free energy G at zero 
temperature in the presence of flow and external pressure 
G = Fo + Ek + PV 

= rrfg + U477g + uqvI H 

with 

^ = ^o + |^ + ^-<^noP + --- 

„, _ _|_ 1 TT^^t** _L 3 m -2 _|_ 

- M04 + 16 + 16 ^-^0 + • • • 

«6 = «06 + i^ + i| + ---, (9) 

where we have expanded Ek also in powers of ?7g and collected the coefficients of 
like powers. Now the details of the supersolid-to-superfluid transition as function 
of the non-ordering parameter, namely, the pressure P and the flow jo> will 
depend on the sign of 1^4 in the present case. (It is to be noted that uq is taken 
to be positive as usual). But quite independently of these details, the effect of 
the zero-point energy and the pressure can be seen clearly from Eq. (9). In 
order to approach the transition, the value of r must decrease sufficiently. This 
is normally prevented here by the largeness of the zero-point term. It is, however. 
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offset by the pressure term so as to bring about the transition. Most significantly 
now, once we are close enough to the transition point, a change in the flow jo can 
drive us across the transition as can be seen from the sign of the term quadratic 
in jo expression for r in Eq. (9) that opposes the pressure term. Indeed, we 
can have then a flow-induced melting of the supersolid. This can be probed 
experimentally. The order parameter rjg determined from the minimization of 
the free energy depends on the flow jo, and enters the non-classical inertial mass 
as in Eq. (8). Inasmuch as for the potential flow in an annular geometry (ID 
flow with periodic boundary condition), the circulation §{jo/no) is quantized to 
vh/m, with v an integer, this inertia will change in corresponding steps as indeed 
observed experimentally [1,2]. 

As for the nature of the transition, in the mean field we have ti4 > give a 
second-order transition. For a first-order transition, we must have < which 
is, of course, unphysical for a helium-like system with inter-atomic repulsive 
interaction. This pathology is due to the ID model that has prevented the 
occurrence of the cubic term in the free energy in Eq. (8). 

Some general remarks are now in order. We have considered here only a 
potential (irrotational) flow and derived the translational non-classical inertia 
associated with the superflow (jo) relative to the density modulation at rest in 
the laboratory frame. (A Galilean transformation can take us to a frame co- 
moving with the modulation relative to the laboratory frame). Experimentally, 
however, it is obviously convenient to have a bounded (confined) motion, which 
is readily realized in rotation — hence the (non-classical) rotational inertia usually 
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measured in experiments by confining the ^He-liquid (solidified under pressure) 
in an annulus which is then made to oscillate about its axis in a torsional mode. 
It is, however, to be noted that for an annular thickness much smaller than the 
annular radius, the motion can be irrotational inasmuch as the annular region 
is not simply connected. Indeed, in the experiments of Kim and Chan [1-2], 
steps corresponding to the Onsager-Feynman quantization of circulation around 
annulus have been seen. Of course, we can have a situation where the motion 
has local circulation distributed in the form of vortices for a simply connected 
system — "^He after all is a type II superfluidi In any case, fundamentally the 
translational inertia is well defined, calculable, and turns out to be non-classical 
as derived above. 

It is apt to recall here that a key point in the microscopic theory of superfluid- 
ity is the role of the condensate (macroscopic occupation of the zero-momentum 
single particle state), namely the hybridization of the single-particle excitation 
and the collective mode caused by the condensate [10]. In the present case, the 
complex modulated order parameter underlying the mass density wave is to be 
viewed as a macroscopically occupied single-particle state, and this too should 
subtend interesting hybridization effects. This calls for further study [11]. 

It may be apt to point out here that the problem of supersolids does raise 
certain general questions of interest about the partitioning of a given amount of 
angular momentum among the different possible modes of motion (degrees of 
freedom) so as to minimize the free energy, or just the energy at zero tempera- 
ture. Thus, e.g., the angular momentum may be shared between the orbital and 
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the spin motion for a system such as in ^He (assuming unpaired spins) giving 
induced spin polarization. In the context of rotating superfluid ^He in a simply 
connected region, it may be taken up by the quantized vortices. Vortices and 
rotons have been invoked recently in the context of supersolids [12]. In a multi- 
ply connected (e.g., annular) region, the angular momentum may betaken up by 
the translational flow as in the present case. For a supersolid spinning about an 
axis, the angular momentum may be taken up by delocalized point defects (the 
defectons where the number of lattice sites exceeds the number of ^He atoms 
[3]. These defectons may have local ring-like exchange motions. For a supersolid 
confined to a thin rotating annulus, however, the quantum defects may carry the 
angular momentum by translating coherently around the annulus. We believe 
that atomistically this may be the case in the recent experiments [1,2]. The 
present single order-parameter based phenomenological theory, however, cannot 
address these atomistic details. 

In conclusion, we have constructed a supersolid model in terms of a com- 
plex order parameter underlying the real mass-density wave characteristic of a 
crystalline solid. This simple model gives a non-classical inertia smaller than the 
literal mass. It also gives the possibility of a flow-induced supersolid-to-superfluid 
transition. Finally, the supersolid as a non-uniform, modulated density n{x) re- 
sults from the fact that underlying the classical looking density modulation there 
is the quantum complex order parameter ijj with = n{x). Thus, a quantum 
supersolid is in a sense the square-root of a classical solid! 

The author would like to thank P. Nozieres for letting him have a copy of his 
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